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ABSTRACT 



m 

Aims. We present a continuation of our study about the relation between stellar mass and gas-phase metallicity in the VIMOS VLT 
^ ■ Deep Survey (VVDS). In this work we extend the determination of metallicities up to redshift x 1.24 for a sample of 42 star-forming 

^—^ • galaxies with a mean redshift value of 0.99. 

' Methods. For a selected sample of emission-line galaxies, we use both diagnostic diagrams and empirical calibrations based on [On] 

, emission lines along with the empirical relation between the intensities of the [Om] and [Nem] emission lines and the theoretical ra- 

tios between Balmer recombination emission lines to identify star-forming galaxies and to derive their metallicities. We derive stellar 
masses by fitting the whole spectral energy distribution with a set of stellar population synthesis models. 

Results. These new methods allow us to extend the mass-metallicity relation to higher redshift. We show that the metallicity de- 
^ ' terminations are consistent with more established strong-line methods. Taken together this allows us to study the evolution of the 

mass-metallicity relation up to z a 1.24 with good control of systematic uncertainties. We find an evolution with redshift of the aver- 
00 ' age metallicity of galaxies very similar to those reported in the literature: for a given stellar mass, galaxies at z ~ 1 have, on average, 

' a metallicity ~ 0.3 dex lower than galaxies in the local universe. However we do not see any significant metallicity evolution between 

redshifts z ~ 0.7 (Paper I) and z ~ 1.0 (this paper). We find also the same flattening of the mass-metallicity relation for the most 
massive galaxies as reported in Paper I at lower redshifts, but again no apparent evolution of the slope is seen between z ~ 0.7 and 
z~ 1.0. 
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1 . Introduction et al. (1979) for dwarf irregular galaxies. Since them, numerous 

works have been performed in order to understand and quantify 



OO 
O 

> 

X 

^ , The formation and evolution of galaxies at different cosmolog- ^^^-^ ^^^^^-^^ ^j^; j^^^l (^mong others SkiUman et al, 
■ - - ■ ical epochs are driven mainly by two processes: the star forma- ^^gg. ^^^^-^ ^ ^^^j^^^^ ^99^. z^ritsky et al., 1994; Richer & 
tion history and the metal enrichment. These processes, which ^^^all, 1994; Garnett et al., 1997 or Pilyugin & Femni, 2000). 
are strongly connected in the evolution of individual galaxies, j^^^ ^^^^^^^^ counterpart in the luminosity-metallicity re- 
can be altered by the pecuhar dynamics of the gas. The rate of ^^^-^^ ^^^-^^^ ^^^-^^ ^^^-^^ ^j^^^ ^ j^^^^ ^^^^^^ ^p^^. 
star formation can be increased as a consequence of mergers or trophotometric bands are available. Besides, in the last years, 
due to the inflow of gas. Regarding metallicity, it is sensitive to ^^^^ ^^^j^ ^^^^ complemented with the use of large surveys 
metal losses due to stellar winds, supernovae and active nucleus g^j^^^^^ ^j^i^j^ ^jj^^ ^j^^ derivation of this relation in sam- 
feedbacks. Therefore, the study of the dominant mechanisms of pj^^ statistically more significant, (e.g. SDSS, Tremonti et al., 
evolution of galaxies at different redshifts rely on the determina- 2004; 2dFGRS, Lamareille et al., 2004), that confirm the de- 
tion of metallicity and total stellar mass. ^^.^^^^ j^^^^j j^^^ ^j^^^^ ^^^^^^ ^^^jj^^. ^^^^ increases. Several 

The correlation between stellar mass and metalhcity (M-Z studies of the evolution of this relation at larger redshifts have 

relation) was firstly found in the nearby universe by Lequeux appeared due to the growing number of spectroscopic surveys 

of star-forming galaxies in the distant universe. This is the case 

Send offprint requests to: E. Perez-Montero of Savaglio et al. (2005) or Lamareille et al. (2006) who agree 
* based on data obtained with the European Southern Observatory fljjd jo^gj. ^ean metallicities for the same galaxy mass bins 
Very Large Telescope, Paranal, Chile, program 070.A-9007, and on j . redshifts. However, they do not agree about the evolu- 

data obtained at the Canada-France-Hawaii Telescope, operated by the ^. r i r^u iv/rr? i t -h » i /^rvrv,r\ 

. _ ^T-TT,^ ■ J J TT ■ K /n ■■ tion of the slope of the M-Z relation: Lamareille et al. (2006) 
CNRS in France, CNRC in Canada and the University ot Hawaii. ,n - ,-i >• i-i ii-,- i 

** n .J . 1 r 11 r »u c ■ u j ■ - - / , reported a flattening of the relation at higher redshifts whereas 

Postdoctoral fellow ot the Spanish Fundacwn Espanola para la , , ^ , .,-,», ^ , • 

Ciencia y la Tecnologfa. Savagho et al. (2005) found a steepening of the M-Z relation. 
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Unfortunately, these studies are limited to z < 1 galaxies, due 
to the required presence of certain emission lines in the optical 
spectrum to derive the metallicity. Erb et al. (2006) studied this 
relation for a sample of star forming galaxies at z ~ 2 using 
near-infrared spectra and they find a lower value of the mean 
metallicity for a given galaxy stellar mass. This evolution is ap- 
parently much more stronger at higher redshift as reported by 
Maiolino et al. (2008), using also IR spectra for a sample of star- 
forming galaxies ar z>3. However, there is a lack of observations 
and available data in the redshift range 0.9 - 2, corresponding to 
the epoch of the universe when the maximum of star formation 
took place. It is therefore particularly interesting to study the M- 
Z relation at this period. 

The estimation of the global trend of the M-Z relation at dif- 
ferent cosmological epochs is possible only through the analysis 
of statistically significant samples with accurate determinations 
of the stellar mass and the gas phase oxygen abundance. This pa- 
per is a continuation of Lamareille et al. (2008; hereinafter Paper 
I), who study the M-Z relation of the star-forming galaxies de- 
tected in the VIMOS VLT Deep Survey up to redshift 0.89. In 
this work we extend the upper redshift limit of the study up to 
z ~ 1 .24, using the relation between metallicity and emission- 
line strength of [On] and [Nem] presented in Perez-Montero et 
al. (2007). In the next section we describe the sample of star- 
forming galaxies discussing the selection criteria and the meth- 
ods to derive their stellar mass and metal content, which can go 
to higher redshift than other classical methods. In section 3, we 
describe the M-Z relation in the redshift bin 0.89 < z < 1.24, 
with an average value in our sample of z = 0.99 and we com- 
pare our results with those found in Paper I. Finally, we sum- 
marise our results and present our conclusions in the last section. 
Throughout this paper we normalize the derived stellar masses 
and the absolute magnitudes with the standard A-CDM cosmol- 
ogy, i.e., h = 0.7, Q,„ = 0.3 and Qa = 0.7 (Spergel et al., 2007). 

2. Description of the sample 

The VIMOS/VLT Deep Survey (VVDS, Le Fevre et al., 2005) is 
one of the widest and deepest spectrophotometric surveys of dis- 
tant galaxies with a mean redshift of z ~ 0.7. The optical spec- 
troscopic data belong to the first epoch sample obtained with 
the Visible MultiObject Spectrograph (VIMOS) at ESO/VLT 
(UT3). They are divided in two deep fields (17.5 < Iab ^ 24): 
VVDS-02h (hereafter F02; Le Fevre et al., 2005) and CDFS (Le 
Fevre et al., 2004) and in three wide fields (17.5 < Iab < 22.5): 
VVDS-lOh, VVDS-14h and VVDS-22h (hereafter F 10, F14 and 
F22 respectively, Garilli et al., in prep.). As explained in Paper 
I, we divide therefore our total sample in two: a deep one with 
the F02 and CDFS fields and a wide one with FIO, F14 and F22 
fields. Spectroscopic data have been reduced using the VIPGI 
pipeline (Scodeggio et al., 2005). 

The following selection criteria have been applied on this 
sample. We consider only objects with a redshift known at a 
75% confidence level {i.e. VVDS redshift flags 2, 3, 4 and 9) and 
we remove all stars, broad-line Active Galactic Nuclei (hereafter 
AGNs; Gavignaud et al., 2006) and duplicate observations. 

The optical spectra of each galaxy have been split into 
two components (stellar and nebular) using the pipeline plate- 
fitjvimos as described in Paper I. This software removes the stel- 
lar continuum and absorption lines and then fits automatically 
all emission lines. The templates used to fit the stellar component 
have been derived from the libraries of Bruzual & Chariot (2003) 
resampled to the velocity dispersion of VVDS spectra. The fit- 
ting and measurement procedures of all emission lines (flux and 
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Fig. 1. VIMOS rest-frame spectrum of the galaxy 
VVDS220485833. The top panel illustrates the fitting of 
the stellar continuum obtained with the platefitjvimos tool. In 
the bottom panel, we show the resultant nebular emission-lines 
spectrum, once the fitted stellar continuum has been substracted. 
The labels indicate the set of emission lines used in this work to 
derive metalUcities. 
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Fig. 2. Histogram of the distribution of the ratio between the 
deredenned fluxes of [Om]4959+5007A and [Neiii]3869A emis- 
sion lines for the five VVDS fields (VVDS-02h, VVDS-lOh, 
VVDS-14h, VVDS-22h and CDFS). The mean value of this dis- 
tribution, shown as a black solid line, is 15.23. The red solid line 
represents the mean value of the distribution without applying 
the mean reddening correction and the dashed line is mean value 
(- 15.37) of this emission-line ratio for a representative sample 
of nearby galaxies studied in Perez-Montero et al. (2007). 



equivalent width) and stellar indices in the original spectra are 
described in Paper I. In Figure [T] we show an example of the fit- 
ting of the stellar continuum (top panel) made by platefit_vimos 
for the galaxy VVDS221485833, one of the best S/N spectra. 
The resultant nebular emission-lines spectrum, once this fitting 
is substracted, is shown in the bottom panel. 
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2.1. Selection of star-forming galaxies 

In order to correctly study the M-Z relation, we select the galax- 
ies whose main mechanism of ionization is star formation. We 
use emission-line diagnostic diagrams to identify narrow-line 
AGNs (Seyfert 2 and LINERS) once broad line AGNs have 
been removed. In particular we use the "blue" diagnostic dia- 
gram based on [Oiii]5007A, Hfi and [On]3727A emission lines 
(Lamareille et al., 2004) and applied in Paper I for the redshift 
bin 0.5 < z < 0.9. We extend this analysis up to redshift ~ 1.3 
using the theoretical ratios between Hfi and the most intense 
Balmer hydrogen recombination lines at bluer wavelengths. This 
is the case of Hy at 4340A and US at 4 102 A. We assume the fol- 
lowing values for these ratios, typical of star-forming galaxies: 
Hr/H/3 = 0.48 and H6/H/3 = 0.25. 

Although these ratios slightly depend on the electron tem- 
perature of the ionized gas, this dependence is negligible in 
comparison with the intrinsic uncertainty of the measurement 
of the involved emission lines, even in the case of narrow-line 
AGNs, which show the lowest values of these ratios. We ap- 
ply the empirical relation found between the emission-line ra- 
tio of [Oiii]4959+5007A and [Neiii]3869A (Perez-Montero et 
al., 2007). The mean value of this ratio reported in this study 
for a sample of HII regions of the local universe is equal to 
15.37. However, the mean value of the same ratio for the VVDS 
star-forming galaxies (both for the deep and wide fields) is 
equal to 17.37. We correct this slight difference taking into 
account a reddening correction corresponding to an extinc- 
tion of Ay = 0.5 mag. After this correction, the mean value 
of [Oiii]4959+5007A/[Neiii]3869A emission-Hne ratio for the 
VVDS galaxies is 15.23, which is almost identical to the value 
found for the reddening corrected HII regions of the local uni- 
verse. Therefore, we decide to use the same reddening correc- 
tion for the rest of the calculations in the VVDS sample. In 
Figure |2] we show the distribution of the emission-line ratio 
[Oiii]4959+5007A/[Neiii]3869A for the 379 VVDS secure and 
candidates star-forming galaxies for which we have measured 
these lines in the redshift range 0.5 < z < 0.9 and with an 
error on the ratio lower than the ratio itself. In this figure we 
compare the mean value of this distribution, once extinction cor- 
rected (black solid line), the mean value found without extinc- 
tion correction (red solid line) and the value found in the local 
universe by Perez-Montero et al. (2007, dashed line). 

The empirical limit between star-forming galaxies and 
narrow-line AGNs in the "blue-Ne" diagnostic diagram based 
on [Nem]3869A emission line takes the following form: 



log 



'[Nell I] 



0.14 



log{[On]IHI3)-l.A5 



•0.23 



(1) 



Following the same prescriptions as in the "blue" diagnostic 
method, the objects lying under this line are considered as star- 
forming galaxies, while those lying above it are considered as 
AGNs. This diagnostic diagram is, in fact, equivalent to that pro- 
posed by Rola et al. (1997) which is based on the same emission 
lines. In the AGN zone, LINERs are identified by the following 
criterium: 



, //([A^e///]3869\ 



-0.57 



(2) 



Adopting the emission-line ratios Hy/Hj6 and WSjWjS listed 
above, we can use the same expression for galaxies at redshifts 
higher than 0.9 or for those whose HyS emission line is not mea- 
sured in the redshift range 0.5 < z < 0.9. In the case of Hd, these 



expressions take the same form as the equations (23) and (24) in 
Perez-Montero et al. (2007). 

In Figure |3] we show the [On]/Hy6 vs. [Neiii]/Hy6 relation for 
the sample of VVDS galaxies with the appropriate lines, in the 
redshift ranges 0.5 < z < 0.9 (left panel) and 0.9 < z < 1.24 
(right panel). The dashed lines delineate the uncertainty domain 
(+0.15 dex around the empirical limit) in the spectral classifica- 
tion: objects lying in this region are considered as "candidates" 
to star-forming galaxies. The results of this analysis are sum- 
marised in Table [T] We make also a comparison with the clas- 
sification derived from the "blue" diagnostic diagram based on 
[On] and [Oiii] emission lines for galaxies in the redshift range 

0. 5 < z < 0.9. 

We see that the fraction of galaxies classified as candidates to 
star-forming from the diagram based on [Nem] is much higher 
(at all redshift bins) that the one based on [Om] emission line. 
Taking into account that the rate of coincidence between blue 
and blue-Ne diagnostics for the galaxies for which it is possi- 
ble to apply at the same time the two methods reaches 98%, we 
think that the sample of galaxies with accurate [Nem] measure- 
ments (i.e. with a high enough S/N) is somehow biased towards 
a population with enhanced [Nem] emission. This is indicated 
by the range of log([Nem]/Hj6) values in both panels of Figure 
[3]which are between -1 and 0, while in Fig. 3 of Paper I we see 
that the ratio log([Om]/HjS), which is approximately an order of 
magnitude larger than log([Nem]/Hy6), varies between -0.5 and 

1. This means that in our sample we are certainly affected by 
some selection effect that preserves us to properly detect all the 
star forming galaxies with a value of log([Nem]/Hy6) between - 
1.5 and -1. Therefore, the relative number of candidates is larger 
in this subsample. 

This same effect is seen when looking at the relative number 
of AGNs in the redshift bin 0.5 < z < 0.9, especially in the 
wide VVDS sample. On the other hand, when we compare the 
relative numbers of AGNs in the redshift regime z > 0.9, we do 
not see a significant increase compared with the blue diagnostic 
diagram at lower redshift. Note, finally, that we did not identify 
any LINER in our VVDS sample. 

For the following analysis of stellar masses and metallici- 
ties, we consider both the samples of candidates and secure star- 
forming galaxies. Unfortunately, there are no deep enough X-ray 
or far-IR observations of the VVDS fields to check the presence 
of hidden AGNs. 



2.2. Estimation of Stellar masses 

The procedure to derive stellar masses is based on the compar- 
ison between observations, basically the photometric Spectral 
Energy Distribution (hereafter SED) of galaxies, combined in 
some cases with their stellar absorption index Hd and D4000 
break measured in VIMOS spectra, and population synthe- 
sis models. As described in Paper I and introduced by SDSS 
collaborators (Kauffmann et al., 2003; Tremonti et al., 2004, 
Brinchmann et al., 2005), we are using the Bayesian approch. 
To summarise, the population synthesis models are based on 
Bruzual & Chariot (2003), with the addition of secondary bursts 
to the standard declining exponential star formation history 
(stochastic library, Salim et al., 2005; Gallazzi et al., 2005). The 
SED of the models are compared with the available photome- 
try for VVDS galaxies: CFH12k observations in BV/?/ bands for 
F02, F22 and FIO fields, and in / band for F14 fields (McCracken 
et al, 2003; Le Fevre et al., 2004), completed by U band observa- 
tions in the F02 field (Radovich et al., 2004), CFHTLS (Canada- 
France-Hawaii Telescope Legacy Survey) in u'g'r'i'z' bands for 
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Table 1. Distribution of star-forming galaxies and narrow-line AGNs among VVDS emission-line galaxies using the "blue" diag- 
nostic diagrams based on the relation between [Oiii] (blue, 0.5 < z < 0.9) or [Nem] (blue-Ne, 0.5 < z < 1.24) and [On] emission 
lines. We show the distribution for the whole sample and for wide and deep fields. In the wide classification we include those 
galaxies belonging to the deep fields having Iab ^ 22.5 
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Fig. 3. Spectral classification of VVDS galaxies using the diagram diagnostic based on [Nem] emission line. In the left panel we 
show the diagram for 418 emission-line galaxies in the redshift range 0.5 < z < 0.9, while the right panel shows 72 emission-line 
galaxies in the redshift range 0.9 < z < 1 .24. In these two panels, the red solid curve divides the regions of star-forming galaxies 
and narrow-line AGNs. The straight red solid line divides the AGN region between Seyfert 2 (top) and LINER (bottom). 



F02 and F22 fields and JKs photometry available for a substan- 
tial fraction of the objects in the F02 field (lovino et al., 2005; 
Temporin et al., 2006). For the CDFS field, we use CFH12k 
UBVRI (Arnouts et al., 2001) and HST bviz (GiavaHsco et al., 
2004) observations. Finally, objects in F02 and F22 fields have 
been cross-matched with the UKIDSS public catolog (Warren et 
al., 2007), providing additional observations in JK bands. 

As in Paper I, we correct for the fact that both minimum and 
maximal stellar detected masses depend on redshift. We thus de- 
fine a limiting mass for each redshift bin, which is the minimum 
mass that would be detected at a given redshift for a given lim- 
iting apparent magnitude. For our M-Z relation study we do not 
consider galaxies with a stellar mass lower that this minimum 
mass derived in the corresponding redshift bin. 



The limiting mass for the redshift range 0.89 < z < 1 .24 is 
derived empirically using the same procedure as in Paper I. We 
derive the limiting mass for each galaxy following the expres- 
sion: 

log(M*.„) = log(M*) + 0.4 ■ (U, - /,,,) (3) 

where M^*^^ is the limiting mass in solar masses, M* the stellar 
mass in solar masses and lohs is the observed /-band magnitude, 
Isei is the /-band magnitude limit for each field, which is equal 
to 24 and 22.5 in the case of deep and wide fields respectively. 
In Figure|4]we show the stellar mass of the selected star-forming 
galaxies as a function of redshift in the range < z < 1.24. 
The solid curve shows the logarithmic fit to the 50% mass-to- 
light ratio completude level of the distribution of the limiting 




redshift redshift 

Fig. 4. Stellar masses of star-forming galaxies (in solar masses) in the VVDS deep (left panel) and wide (right panel) fields as a 
function of redshift. The squares represent the star-forming galaxies selected using the diagnostic diagram based on [Neiii] emission 
line. The solid curve represents a logarythmic fit to the 50% mass-to-light completude levels and shows the bias imposed by the 
VVDS selection on apparent magnitude. 
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Fig. 5. Comparison between O/H abundances derived for the 
star-forming galaxies of the SDSS DR4 catalogue using the 
02A'f3 parameter under different assumptions about reddening. 
Filled histogram represents the difference between the abun- 
dance using emission lines corrected for reddening using the 
Balmer decrement and assuming an average reddening corre- 
sponding to Ay - 0.5 mag. The thick line histogram represents 
the same difference but in relation with the abundances derived 
from the equivalent widths. We indicate also the mean value and 
the standard deviation of the two distributions. 

mass in each redshift bin. Since the average redshift value for 
the analysed galaxies in the range 0.89 < z < 1.24 is z = 1.0, we 
take a value for the limiting mass in this regime of log(M»/M0) 
= 8.9 in the deep fields and 9.6 in the wide fields. As in Paper I, 
we obtain a lower value of the limiting mass in the deep fields. 
The values are also slightly larger than those derived for the star- 




Fig. 6. Relation between the 02Ne3 parameter and the emission- 
line ratio [NeIII]/[OII] for the VVDS star-forming galaxies at 
0.9 < z < 1.24. OiNei has been calculated using either the flux 
of Hj3 (blue squares), Hy (black squares), or H6 (white squares). 



forming galaxies in the redshift bin 0.7 < z < 0.9 which are 8.8 
and 9.5 respectively. This is due to the fact that limiting mass 
increases with redshift in the selection function. 

2.3. Estimation of metallicities 

We estimate metallicities as the oxygen abundances of the ion- 
ized gas phase from the 02Ne3 parameter (Perez-Montero et al., 
2007). This parameter is equivalent to O23 (also known as R23), 
but it is based on the empirical relation existing between [Oiii] 
and [Neiii] emission lines: 

/([0//]3727A) H- 15.37 ■ /([A^e///]3869A) ^ 
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12+log (O/H) (02Ne3 (HP)) 12+log (O/H) (02Ne3 (Hy)) 12+log (O/H) (02Ne3 (H5)) 

Fig. 7. Comparison between the oxygen abundances derived with the O23 and 02Nei parameters. Three cases are distinguished: 
02Ne3 parameter is estimated using H/3 (left panel, 268 galaxies), Hy (middle panel, 217 galaxies) or Hd (right panel, 187 galaxies) 
emission line. Only VVDS star-forming galaxies with an accuracy on O/H lower than 0.3 dex are considered. Dashed lines represent 
the 1 : 1 relation while the solid lines are the best linear fit to the data. 



This calibration allows to derive metallicities using the theo- 
retical ratios between H/3 and the other bright hydrogen recombi- 
nation lines at bluer wavelengths, like Hy and H6. It is therefore 
possible to derive metallicities of galaxies up to redshift z ~ 1 .25 
with the optical VIMOS spectra obtained by the VVDS. In that 
case, the 02Nei parameter takes the following expression: 

_ /([(9//])+ 15.37 -/([jVg///]) KHi) 

^^'^^^ im IW) ^ 

We keep the factor 15.37 obtained in Perez-Montero et al. 
(2007) for a sample of HII regions of the local universe once the 
observed fluxes have been dereddened. As shown in Figure |2] 
we obtain a very similar factor to correct the flux of the involved 
emission lines, with a reddening equivalent to an extinction of 
Ay=0.5 mag. This correction has a very low impact on the de- 
termination of metallicities obtained with this method: they are 
only 0.03 dex lower for those derived using Hy and 0.01 dex 
lower using Hd. Contrary to Paper I where equivalent widths 
were used as proposed by Kobulnicky & Phillips (2003) and 
Lamareille et al. (2006a), emission line fluxes work better in 
the case of 02Nei parameter because the variations of the con- 
tinuum are more sensible to the underlying stellar population in 
this part of the spectrum. To illustrate this point we have done 
an analysis of the star-forming galaxies of the DR4 SDSS cat- 
alogue (Kauff'man et al., 2003) whose involved emission lines 
have a signal-to-noise larger than 5. We show in Figure |5] the 
histograms of the residuals to the oxygen abundances derived 
from 02Ne3 and whose emission lines have been corrected for 
reddening using the constants derived from the B aimer decre- 
ment. The filled histogram represents the residuals in relation 
with the abundances using the same parameter, but with an av- 
erage reddening correction for all the sample equivalent to Ay 
= 0.5 mag. On the other hand, the empty thick line histogram 
represents the residuals with the abundances deduced with the 
same parameter, but using equivalent widths and correction fac- 
tors based on the shape of the continuum calculated in a similar 
way as in Kobulnicky & Phillips (2003). As we can see, both the 
mean value of the residuals and the dispersion of the distribution 
are much better in the first case. 

As described in Perez-Montero et al. (2007), the 02Aff3 pa- 
rameter has the same problems than O23, the most important be- 
ing its degeneracy with metallicity. This means that for a given 
value of the 02Ne3 parameter we can find two possible values of 
the metallicity. Among the numerous calibrations of the O23 pa- 



rameter there is always a calibration for high metallicities (upper 
branch) and another one for low metallicities (lower branch) (e.g. 
Skillman, 1989). It is thus mandatory to be able to distinguish 
between these two metallicity regimes. The ratio [Neiii]/[On] 
has been proposed by Nagao et al. (2006) as an empirical cal- 
ibrator of the metallicity but, due to its high dependence on ion- 
ization parameter, it does not allow to retrieve accurate deter- 
minations of O/H. However, it can be used to choose the ap- 
propriate branch (low or high metallicity regime) for the 02Ne3 
parameter The results obtained by Perez-Montero et al. (2007) 
indicate that galaxies with log([Neiii]/[Oii]) < -1.0 are in the 
high metallicity regime, although the contrary is uncertain. In 
Figure |6] we show the relation between 02Ne3 parameter and 
log([Nem]/[Oii]) for the VVDS galaxies harbouring the involved 
emission lines. We see that only a fraction of the objects lie in 
the high metallicity regime using this criterion. For the rest of 
galaxies with log([Neiii]/[Oii]) > -1.0, we constrain the metal- 
licity regime by comparing the derived metallicities in low- and 
high-metallicity branches to the value derived with the Bayesian 
method described in Paper I or with the O23 parameter in the 
case of the redshift bin 0.5 < z < 0.9. After this comparison we 
conclude that 95% of the sample of star-forming galaxies lie in 
the high metallicity regime. The selected objects have values of 
02Ne3 in the range 0.0-1.5, corresponding to oxygen abundances 
in the range 8.0 < 12-Hlog(0/H) < 9.2 in the upper branch. 

In order to be consistent with metallicity estimates per- 
formed in Paper I, we use the calibrations of the upper branch 
of O23 parameter based on the models of Chariot & Longhetti 
(2001). As it has been shown by Kewley & Eflison (2008), it 
is very important to keep the consistency between the different 
empirical calibrators used to derive the metallicity of galaxies. If 
not, systematic differences as high as 0.7 dex can be observed for 
the M-Z relation in the local universe as derived from the SDSS. 

We check the consistency of this new method at higher red- 
shift by comparing the metallicities derived using the 02Ne3 and 
O23 parameters in the VVDS star-forming galaxies at z < 0.89. 
Up to this redshift, all the emission lines required to compute 
these two parameters are observed in the VIMOS optical spec- 
tra. For this purpose we take into account only the oxygen abun- 
dances derived with an accuracy better than 0.3 dex. The re- 
sults of this comparison are shown in Figure |7] We distinguish 
three cases where 02Ne3 parameter is estimated using Hfi (268 
galaxies), Hy (217 galaxies) and H6 (187 galaxies) respectively. 
Dashed lines represent the 1 : 1 relation while the solid lines are 
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the best linear fit to the data. The best correlation is found in the 
first case, where the average of the residuals is 0.01 dex with 
a standard deviation of 0.17 dex. For the second case (middle 
panel) the mean value of the residuals is also 0.01 dex, but the 
standard deviation increases to 0.25 dex. Finally, in the third case 
(right panel), we obtain a mean value of the residuals of 0.04 dex 
and a standard deviation of 0.24 dex. The same kind of analysis 
comparing the metallicities derived from the Oawes parameter on 
one side and the Bayesian Method on the other side gives a mean 
value of the residuals of 0.03 dex only with a standard deviation 
of 0.31 dex. 

3. The mass-metallicity relation 

We study the M-Z relation for the redshift bin 0.89 < z < 1.24, 
using the stellar masses and the metallicities derived in previ- 
ous sections. We take only galaxies with a stellar mass higher 
than the limiting mass imposed in Figure |4] for a mean red- 
shift z = 1, corresponding to a logarithmic stellar mass in solar 
masses of 8.9 for the deep fields and 9.6 for the wide fields. We 
take only galaxies with an uncertainty in the calculation of the 
oxygen abundance lower than 0.3 dex, so they remain 42 star- 
forming VVDS galaxies: 27 in the deep fields and 15 in the wide 
fields. The galaxies are listed in Table|2]along with their Iab ap- 
parent magnitude, spectroscopic redshift, and derived metallic- 
ities and stellar masses. In the case of metallicity, all the errors 
due to the uncertainties on the emission-line fluxes are lower 
than the intrinsic uncertainty of the empirical parameter based 
on [Neiii] emission lines, which is 0.2 dex for the calibration 
of the upper branch of the 02Ne3 parameter (Perez-Montero et 
al., 2007). The spectra of all these galaxies are available via the 
CENCOS databas43- We have to notice that, after applying the 
selection criteria described above, no galaxies from the lOh and 
14h VVDS fields remain in the final sample. 

In left panel of Figure[8]we show the relation between stellar 
mass and oxygen abundance for these 42 galaxies. For compar- 
ison, we show also the M-Z relation derived in the local uni- 
verse (z < 0.3) by Tremonti et al. (2004) with SDSS data (solid 
line). Note that the comparison can be made directly as the same 
method/calibration (namely Chariot & Longhetti 2001, see sect. 
2.3) has been used to derive metallicities in the local and high- 
redshift galaxies. The respective limiting mass for both deep and 
wide fields are shown as vertical solid lines. 

In order to test the global evolution in metallicity of star- 
forming galaxies we make the same analysis carried out in Paper 
I, consisting in the calculation of the mean shift on metallicity by 
fitting the same curve deduced for the SDSS sample to our data 
points and, hence, the derivation of a different zero-point. These 
fits are shown in left panel of Figure [8] as a blue dashed-point 
line for the galaxies in the deep fields and a red dashed-point 
line for the wide field galaxies. The mean difference in metallic- 
ity found with respect to the local M-Z relation are -0.29+0.21 
and -0. 12+0. 15 for the wide and deep field galaxies respectively. 
As it was shown in Paper I, we see a stronger evolution of the 
mean metallicity of galaxies in the wide sample. As the limiting 
mass of the wide sample is much higher than in the deep sample, 
this could be due to a a more significant evolution in metallic- 
ity for the most massive galaxies. The M-Z relation derived in 
the redshift range 0.89 < z < 1.24 is rather similar to the one 
derived in Paper I for the redshift bin 0.7 < z < 0.9. However, 
this difference between the wide and deep samples could be due 
to different selection effects. First the exposure time was higher 

' |http : // cencosw . oamp . f r | 



for the deep fields which could lead to a higher dispersion in the 
metallicity distribution of these objects, as the number of high- 
metallicity galaxies, i.e. with fainter emission lines, is larger 
There is also a larger number of low-mass galaxies in the deep 
sample, for which the metallicity is also lower compared with 
the average value. 

We assess also the shape of the M-Z relation calculating a 
running median of the metallicity for the whole sample (both 
deep and wide fields) along the stellar mass axis in order to com- 
pare it with the results from Paper I and the relation deduced 
for the local universe using the SDSS sample (Tremonti et al. 
2004). The method to calculate these median is identical to that 
described in Paper I: we have taken four bins in stellar mass, 
with a similar number of objects in each bin («; 1 1) and then we 
calculate the median of 1000 bootstraps estimates. These medi- 
ans are shown in the right panel of Figure [8] as black squares. 
The error bars associated to the median metallicities of galaxies 
with z > 0.9 are much higher than the ones derived in Paper I at 
lower redshifts, due to the lower number of objects in each mass 
bin. The M-Z relation found for 0.9 < z < 1.24 is fully com- 
patible with the results of Paper I for z > 0.5 in absolute value 
and shape, pointing to a significant global chemical evolution for 
redshifts between and 0.5. However, for a redshift z > 0.5, we 
do not find a significant global evolution within the error bars 
up to z = 1 . Regarding the shape of the M-Z relation, we see 
a flattening as compared the with the M-Z relation in the local 
universe. This could be a consequence of more effective yields 
for the more massive galaxies, where the evolution is stronger. 
However, again, we do not appreciate any significant evolution 
of the slope of the M-Z relation between redshifts 0.5 and 1. 

It is expected somehow that the sample of star-forming 
galaxies considered in this work can be affected by some selec- 
tion effects which lead to a higher fraction of objects with low 
metallicities. Firstly, the selection criterion on [Nein] emission 
lines biases the sample towards lower metallicities as pointed by 
Nagao et al. (2006) and Perez-Montero et al. (2007). Secondly, 
the criterion on hydrogen recombination lines (Hy and H6) im- 
plies a sample of galaxies with higher star-formation efficiencies. 
According to Ellison et al. (2008), this leads also to galaxies with 
lower metallicities. Nevertheless, the consistency between the 
metallicities derived using the 02Ne3 parameter in this redshift 
range and those at lower redshifts indicates that the probable se- 
lection effects inherent to our sample do not have a significant 
weight in the overall deduced metallicity. 



4. Summary and conclusions 

We have carried out a continuation of our study of the mass- 
metallicity relation of star-forming VVDS galaxies extending 
our analysis up to redshift z ~ 1 .24. 

We selected galaxies with a reliable measurement of [On], 
[Neiii] and Hy or H6 emission lines and we identified star- 
forming galaxies with the diagnostic diagram based on these 
lines and presented in Perez-Montero et al. (2007). Then, we de- 
rived oxygen gas-phase abundances using the empirical parame- 
ter OiNei- We have shown that these new methods are consistent 
with those based on [Oiii] if we use the same calibrations. As it is 
based on weaker emission lines that, generally, are more affected 
by absorption effects, the only problem appears as a selection ef- 
fect, biaising the sample towards galaxies with bright emission 
lines and, hence, towards lower metallicities. The relative num- 
ber of narrow-line AGNs at higher redshifts is thus to be taken 
with caution. 
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Table 2. Sample of 42 star- forming VVDS galaxies in the redshift regime 0.89 < z < 1.24 with a determination of both the stellar 
mass and the oxygen abundance. We report also their spectroscopic redshift and /-band (AB system) apparent magnitude. 



Object 




z 


12+log(0/H) 


log(M./Mo) 


Object 




z 


12+log(0/H) 


log(M./Mo) 


VVDS000016597 


23 


24 


0.9689 


8.68 


9.24±0.30 


VVDS020386743 


22.58 


1.0465 


8.91 


9.43+0.11 


VVDS000025365 


22 


80 


0.8956 


8.50 


10.10±0.25 


VVDS020389342 


23.76 


0.9057 


8.36 


9.30+0.14 


VVDS000028748 


21 


24 


0.9523 


9.16 


9.42±0.29 


VVDS020460957 


23.36 


0.9615 


8.96 


9.34+0.08 


VVDS020098524 


23 


13 


0.9044 


8.68 


8.95+0.07 


VVDS020461893 


23.45 


1.0486 


8.73 


9.27+0.07 


VVDS020109569 


22 


50 


1.0800 


8.51 


9.40+0.05 


VVDS020462475 


22.76 


0.9626 


8.35 


9.36+0.07 


VVDS020 127928 


21 


23 


0.9610 


8.89 


10.43+0.07 


VVDS220005349 


22.36 


0.9610 


9.03 


10.20+0.14 


VVDS020 142000 


22 


16 


0.9606 


8.52 


9.45+0.04 


VVDS220077555 


22.10 


0.9825 


8.25 


9.89+0.07 


VVDS020 150832 


23 


19 


0.9628 


8.38 


9.16+0.06 


VVDS220102663 


21.83 


1.1348 


9.01 


10.14+0.11 


VVDS020155109 


22 


17 


0.9840 


8.57 


9.62+0.06 


VVDS220108994 


22.3 


0.9007 


8.16 


9.36+0.11 


VVDS020 165728 


23 


11 


0.9051 


8.62 


9.06+0.05 


VVDS220156608 


22.06 


0.9871 


8.68 


10.53+0.11 


VVDS020175617 


23 


18 


0.9046 


8.72 


8.94+0.16 


VVDS220294849 


22.17 


0.9933 


8.65 


9.93+0.17 


VVDS020 196604 


23 


12 


0.9607 


8.60 


9.06+0.14 


VVDS220376206 


21.85 


1.2440 


8.74 


10.02+0.08 


VVDS020199214 


22 


62 


0.9126 


8.27 


9.28+0.07 


VVDS220380153 


21.55 


0.9332 


8.72 


9.94+0.06 


VVDS020201361 


23 


27 


1.0572 


8.83 


9.03+0.14 


VVDS220390996 


21.76 


0.9071 


8.41 


9.75+0.09 


VVDS0202 10745 


23 


97 


0.8983 


8.61 


9.31+0.14 


VVDS220457904 


22.30 


1.0666 


8.48 


10.06+0.08 


VVDS020247631 


23 


37 


0.9171 


9.04 


9.03+0.10 


VVDS220485833 


20.99 


0.9807 


8.58 


10.33+0.10 


VVDS020276040 


21 


71 


0.9249 


9.02 


10.05+0.07 


VVDS220535047 


22.29 


1.0159 


8.67 


10.01+0.10 


VVDS020294045 


22 


80 


1.0019 


8.29 


9.17+0.06 


VVDS220540780 


21.91 


0.9089 


8.72 


10.85+0.14 


VVDS020304000 


23 


46 


0.9561 


8.64 


9.28+0.13 


VVDS220559147 


22.25 


0.8984 


8.99 


10.48+0.14 


VVDS020322281 


22 


89 


0.9320 


8.65 


9.32+0.07 


VVDS220603467 


21.85 


0.9156 


8.81 


10.31+0.08 


VVDS020385905 


22 


36 


0.9759 


8.44 


9.84+0.01 


VVDS220605155 


22.06 


0.9519 


8.52 


9.78+0.07 




8,0 



J_ 



• deep fields 
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Fig. 8. Relation between the gas-phase metallicity and the stellar mass (in solar masses) for VVDS star-forming galaxies in the 
redshift range 0.89 < z < 1.24. In left panel, we show the 42 individual VVDS star-forming galaxies (black circles: deep fields; 
white circles: wide fields). Vertical solid lines represent the corresponding limiting mass to the deep and wide samples. The solid 
red line represents the Tremonti et al. (2004) M-Z relation for the SDSS sample. We show as well the typical uncertainty of the 
02^63 parameter as vertical eiTor bars. The blue dashed-point curve represents the fit of the Tremonti curve to the VVDS galaxies of 
the deep sample, and the red for the wide sample. In right panel, we show the running medians of the whole sample (black squares, 
both deep and wide fields), as compared with the results from Paper I at lower redshifts (coloured squares). 



The final sample used to derive the relation between metal- 
licity and stellar mass consists of 42 VVDS star-forming galax- 
ies with uncertainties lower than 0.3 dex in both oxygen abun- 
dance and stellar mass. These have been selected to be larger 
than the previously calculated limiting mass for this redshift 
range. Since the number of selected objects is sensibly lower 
than in our previous study (Paper I), the sample is not as statisti- 
cally significant but we have found fully consistent results. 



The zero-point and the slope of the M-Z relation in the 
0.9 < z < 1 .24 redshift range for both the deep and wide fields 
are rather similar to the values reported in Paper I for z > 0.5, 
showing a significant global cosmic chemical evolution in rela- 
tion with the local universe. We detect also a flattening of the 
M-Z relation in comparison with the SDSS-based relation in the 
local universe (Tremonti et al., 2004). This evolution is stronger 
for the galaxies of the wide sample, consistently with the fact 
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that most massive galaxies have undergone a stronger chemical 
evolution. The evolution of massive galaxies compared with the 
local sample can be understood as a consequence of lower ef- 
fective yields for the most massive galaxies at larger redshifts, 
according to an hierarchical scenario of formation and evolution 
of galaxies. 

Nevertheless, despite the probable selection effects of the fi- 
nal sample, we do not appreciate any significant evolution nei- 
ther in the absolute value of the metallicity nor in the shape 
of the M-Z relation for redshifts larger than 0.5. This can be a 
consequence of the smaller time interval in the redshift range 
0.7 < z < 1 combined with a larger uncertainty due to the lower 
number of objects with an appropriate signal-to-noise ratio for 
the used emission-lines. In any case, our results are fully com- 
patible with those published for z ~ I in the DEEP2 sample 
(Shapley et al., 2005; Liu et al., 2008) or in the GOODS sample 
(Cowie & Barger, 2008), who find 0.2-0.3 dex lower metallici- 
ties in relation with the local universe. 
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